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The variable magnetic field of the solar photosphere exhibits periodic reversals as a result of dynamo activity occur- 
ring within the solar interior We decompose the surface field as observed by both the Wilcox Solar Observatory and the 
Michelson Doppler Imager into its harmonic constituents, and present the time evolution of the mode coefficients for the 
past three sunspot cycles. The interplay between the various modes is then interpreted from the perspective of general 
dynamo theory, where the coupling between the primary and secondary families of modes is found to correlate with 
large-scale polarity reversals for many examples of cyclic dynamos. Mean-field dynamos based on the solar parameter 
regime are then used to explore how such couplings may result in the various long-term trends in the surface magnetic 
field observed to occur in the solar case. 

Subject headings: Sun: interior — Sun: photosphere — Sun: magnetic fields — Sun: activity 



1. Introduction 

The Sun is a dynamic star that possesses quasi-regular cy- 
cles of magnetic activity having a mean period of about 22 yr. 
This period varies from cycle to cycle, and over the past sev- 
eral cen turies has ranged from 18 to 25 years (IWeiss 1 19901: 
Beer et al. 1998; Usoskin et al. 2007J, as for example illus- 
trated by the unusual but not unprecedented length of the most 
recently completed sunspot cycle 23. During each sunspot 
cycle (comprising half of a magnetic cycle), the Sun emerges 
sunspot groups and active regions onto the photosphere, with 
such features possessing characteristic latitudes, polarity, and 
tilt angles. As with the period, the numbers and emergence 
frequencies of active regions is observed to vary from cycle 
to cycle. 

At activity minima when few active regions are present, 
the surface magnetic field is characterized by the presence 
of two polar caps, i.e., largely unipolar patches of magnetic 
flux dispersed across both polar regions with the northern and 
southern caps possessing opposite polarities. Reversals of this 
large-scale dipole represented by the polar-cap flux occur dur- 
ing each sunspot cycle, allowing the subsequent sunspot cy- 
cle to begin in the opposite configuration. After two sunspot 
cycles, and thus after undergoing two polarity reversals, the 
photospheric field will have returned to its starting configura- 
tion so as to complete a full activity cycle. 

In response to the photospheric flux associated with vari- 
ous features, such as active regions and their decay products, 
the coronal magnetic field possesses structures having a broad 
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spectrum of sizes. These structures are both evident in obser- 
vations of coronal loops, as found in narrow -band extreme ul- 
traviolet or soft X-ray imagery, a nd reproduced in models o f 
the coronal magnetic field (e.g., ISchriiver & DeRosa 2003 ). 
In both venues, the coronal magnetic field geometry is seen 
to contain a rich and complex geometry. Dynamical events 
originating from the corona, such as eruptive flares and coro- 
nal mass ejections, are likely powered by energy released by 
a reconfiguration of the coronal magnetic field, which in turn 
is responding to changes and evolution of photospheric fields. 

Precise measurements of the time-history of photospheric 
magnetic field and the ability to determine the projection 
of this field into its constituent multipole components are 
helpful in investigating the physicaljjrocesses thoiight to be 
responsible fo r dynamo activity feullard & GellmanI 1 195i 
MoffattI 1978h . In cool stars similar to the Sun, the dynamo 
is presumed to be a consequence of the nonlinear interactions 
between convection, rotation, and large scale flows, lead- 
ing to the generation and maintenance against Ohmic diffu- 
sion of magnetic field of various temporal and sp ati al scales 
(IWeisslll987l: ICattaneolll999l: IOssendriiveill2003l: iBnin et al 



200 4 IVogler & Schussleill2007t ICharbonneaull2010l: iReiners 
2012h . In particular, the dependence of dy namo activity upon 
rotation appears to be well estabhshed (Novesetal. 1984t 



Saar & Brandenburg' 1999 HPizzolato et all2003l:lBohm-Vitense 
2007; Reiners et al. 2009|). However, many details of the un- 
derstanding of why many cool-star dynamos excite waves of 
dynamo activity having a regular period, specifically 22 yr in 
the case of the Sun, remain unclear. 

To investigate this question, it is useful to explore the be- 
havior and evolution of the lowest-degree (i.e., largest-scale) 
multipoles, their amplitudes and phases, and thek correla- 
tions with the solar photospheric magnetic field. Many earher 
studies (e.g., lLevinelll977l ; iHoeksemal 1 1 9841 ; iGokhale et al 
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1992l:lGokhale & Javaraiahll 19921) have illustrated how power 
in these modes ebb and flow as a function of the activity level. 
In particular, J. Stenflo and collaborators have performed 
thorough spectral analyses on the t emporal evolution of th e 



various spherical harmonic m odes. Stenflo & Voge| 
Stenflo & WeisenhornI (fl987l) and Ist enflo & Giidel 



19861 



19881 



and m ore recently Stenflol ( 1994 ) and Knaack & Stenflo 



(120051) . base their analysis on Mt. Wilson and Kitt Peak mag- 
netic data spanning the past few sunspot cycles. As one would 
expect, they find that most of the power is contained in tem- 
poral modes having a period of about 22 yr, and especially 
in spherical harmonics that are equatorially antisymmetric, 
such as the axial dipole and octupole. However, they find 
signatures of the activity cycle are present in all axisymmet- 
ric harmonics, as significant power is present at temporal 
frequencies at or near integer multiples of the fundamental 
frequency of 1.44 nHz [equivalent to (22 yr)^^]. 

In the current study, we focus on the coupling between 
spherical harmonic modes, and what such coupling may indi- 
cate about the operation of the interior dynamo. In particular, 
reversals of the axial dipole mode may be viewed as a result 
of continuous interactions between the poloidal and toroidal 
components of the interior magnetic field, i.e., the so-called 
dynamo loop. Currently, one type of solar dynamo model 
that successfully reproduces many observed behaviors is the 



flux-t r ansport Babcock-Leighton tvpe (e g.. iChoudhuri et al 



1995 



Dikpatietal.1 120041: iJouve & BrunI 120071: lYeates et al 



2008r) . A key ingredient in producing realistic activity cycles 
using this type of model is found to b e the amplitude and pro- 
file of the meridional flow (Ijouve & B run 2007; Karak 2010t 
Nandy et al.ll201 ll: lDikpatill2oT Ij), which result in field rever- 
sals progressing via the poleward advection across the sur- 
face of trailing-polarity flux from emergent bipolar regions. 
During the rising phase of each sunspot cycle, polar cap flux 
left over from the previous cycle is canceled, after which 
new polar caps haying the opposite magnetic polarity form 



( Wang et al. 19891 : Benevolenskayal 2004 : Dasi-Espuig et al 
2010h . 



Helioseismic analyses of solar oscillations have provided 
measurements and inferences of key dynamo components, 
such as the internal rotation profile and the near-surface 
meridional circulation (Thompson et al. 2003; Basu & Antia 



20101) . Complementing precise observations of the solar mag- 
netic cycle properties, these helioseismic inversions repre- 
sent additional strong constraints on theoretical solar dynamo 
models. Successful solar dynamo models strive to reproduce 
as many empirical features of solar magnetic activity as pos- 
sible, including not only cycle periods, but also parity, phase 
relation between poloidal and toroidal components, and the 
phase relation between the dipole and higher-degree harmonic 
modes. 

Interestingly, a recent analysis of geomagnetic records has 
indicated that the interplay between low-degree harmonic 
modes during polarity reversals is one way to characterize 
both reversals of the geomagnetic dynamo (which have a 



mean period of about 300,000 yr) as well as excursions, 
where the dipole axis temporarily moves equatorward and 
thus away from its usual position of being approximately 
aligned with the rotation axis, followed by a return to its 
original position without having crossed the equator (see 
Hulot etaLlEoiol for a recent review on Earth's magnetic 



field). In particular, these studies have shown that, during pe- 
riods of geomagnetic reversals, the quadrupolar component 
of the geomagnetic field is stronger than the dipolar compo- 
nent, while during an excursion (which can be th ought of as 
a failed reversal), the dipole remains dominant ( Amit et al 



20101: iLeonhardt & Fabianll2007tlLeonhardt et alj|2009h . One 

may thus ask: Is a similar behavior observed for the solar 
magnetic field? 

In an attempt to address this question, we have performed 
a systematic study of the temporal evolution of the solar pho- 
tospheric field by determining the spherical harmonic coeffi- 
cients for the photospheric magnetic field throughout the past 
three sunspot cycles, focusing on low-degree modes and the 
relative amplitude of dipolar and quadrup olar compon ents. 
Following the classification of McFadden et al.l ( 1991 ). we 
have made the distinction between primary and secondary 
families of harmonic modes, a classification scheme that takes 
into account the sym metry and parity of the spherical har- 
monic functions (see iGubbins & Zhanglll993l for a detailed 
discussion on symmetry and dynamo). 

While we recognize that the solar dynamo operates in a 
more turbulent parameter regime than the geodynamo, and is 
more regular in its reversals, the presence of grand minima 
(such as the Maunder Minimum) in the historical record in- 
dicates that the solar dynamo can switch to a more intermit- 
tent state on longer-term, secular time scales. In fact in the 
late stages of the Maunder Minimum, the solar dynamo was 
apparently asymmetric, with the sou thern hemisphere pos- 
sessing more activity than the north jRibes & Nesme-Ribes 
[l993 ) for several decades, a magnetic configuration that may 
have been achieved by having dipolar and quadrupolar mode s 
of similar amplitude ("Tobias! 1 1997t ICaUet & Petrelisll2009h . 
Additionally, recent spectopolarimetric observations of solar- 
like stars now provide sufficient resolution to characterize the 
magnetic field geometry in terms of its multipolar decom- 
position (Petit et al. 2008). Furthermore, the analysis of re- 
duced dynamical systems developed over the last 20 yr de- 
scribing the geodynamo and solar dynamo have emphasized 
the importance of the nonlinear coupling between dipolar 
and quadrupolar components iKnobloch & Landsberglll996t 



Weiss & Tobiadl2000t IPetrehs et al.11200 9'). 

This article is organized in the following manner In ^ 
we describe the data sets and the data analysis methods used 
to perform the spherical harmonic analysis, followed in ^ 
with an explanation of the temporal evolution of the various 
harmonic modes, the magnetic energy spectra, and the de- 
composition in terms of primary and secondary families. We 
interpret in ^our results from a dynamical systems perspec- 
tive and illustrate some of these concepts using mean-field 
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dynamo models. Concluding remarks are presented in ^ 

2. Observations and Data Processing 

We analyze time series of synoptic photospheric mag- 
netic field maps of the radial magnetic field Br derived from 
line-of-sight magnetogram observations taken by both the 
Wilcox Solar Observatory (WSO; Scherrer et al. 197^ at 
Stanfo rd University and by the Michelson Doppler Imager 
(MDI: IScherrer et al.lll995h on board the space-borne Solar 
and Heliospheric Observatory (SOHO). The WSO datiflused 
in this study span the past 36 years, commencing with Car- 
rington Rotation (CR) 1642 (which began on 1976 May 27) 
and ending with CR 2123 (which ended on 2012 May 25). 
For MDo we used data from much of its mission lifetime, 
starting with CR 1910 (which began on 1996 Jul. 1) through 
CR 2104 (which ended on 2010 Dec. 24). In both data series, 
one map per Carrington rotation was used, though maps with 
significant amounts of missing data were excluded. The mea- 
sured line-of-sight component of the field is assumed to be 
the consequence of a purely radial magnetic field when calcu- 
lating the harmonic coefficients. Additionally, for WSO, the 
synoptic map data are known to be a factor of about 1.8 too 
low d ue to the saturation of the instrument (iSvalgaard et al 



19781) . Lastly, the MDI data have had corrections applied for 
the polar fields using the interpolation scheme presented in 
SunetalJ(l201lb . 

For each map, we perform the harmonic analysis using 
the Legendre-transform software provided by the "PFSS" 
package available through SolarSoft. Using this software 
first entails remapping the latitudinal dimension of the in- 
put data from the sine-latitude format provided by the ob- 
servatories onto a Gauss-Le gendre grid (c.f., § 25.4.29 of 
Abramowitz & Stegun Il972h . This regridding enables Gaus- 
sian quadrature to be used when evaluating the sums needed to 
project the magnetic maps onto the spherical harmonic func- 
tions. The end result is a time-varying set of complex coeffi- 
cients (t) for a series of modes spanning harmonic degrees 
£ = 0, 1, £rnax, where the truncation limit ^max is equal to 
60 for the WSO maps and 192 for MDI maps. The co- 
efficients are proportional to the amplitude of each spherical 
harmonic mode Y^™ for degree £ and order m possessed by 
the time series of synoptic maps, so that 



(1) 



where 9 is the colatitude, cf) is the latitude, and t is time. We 
note that because the coefficients are complex numbers, 
this naturally accounts for the rotational symmetry between 
spherical harmonic modes with orders m and — m (for a given 
value of £), with the amplitudes of modes for which to > 
appearing in the real part of i?™^ ™d the amplitudes of the 



-'Available at'ht tp :/ / wso ■ Stanford ■ edu / synopticl ■ html' 
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modes where to < being contained in the imaginary part 
of B"\ Consequently, the sum over ni in equation ([T]i starts 
at TO = instead of at m = The coefficients corre- 
sponding to the axisymmetric modes (for which to = 0) are 
real for all £. 

We use the convention that, for a particular spherical har- 
monic degree £ and order to. 



(2) 



where the functions P™(cos(?) are the associated Legendre 
polynomials, and where the coefficients C™ are defined 



CI' 



i-iy 



2£+l {£-m)\ 
Att {£ + m)l 



(3) 



With this normalization, the spherical harmonic functions sat- 
isfy the orthogonality relationship 



sin 9 dOrr 



= 5u' 5, 



(4) 



When comparing our coefficients with those from other stud- 
ies, it is important to take the normalization into account. For 
example, the complex coefficients used here are different 
than (albeit related to) the real-valued 5™ and /i™ coefficients 
provided by the WSO tearr0. This difference is due to their 
use of spherical harmonics having the Schmidt normalization, 
a convention that is commonly used by the geomagnetic com- 
munity as well as by earlier studies in the solar community 
such as Altschuler & Ne wkirk ( 1969). For the WSO data used 
here, we have verified that the values of used in this study 
are commensurate with the and ft,™ coefficients provided 
by the WSO team. 

Because we possess perfect knowledge of Br neither over 
the entire Sun nor at one instant in time, the monopole co- 
efficient function B'^{t) does not strictly vanish and instead 
fluctuates around zero. In practice, we find that the magni- 
tude of B^{t) is small, and thus feel justified in not consid- 
ering it further. This assumption effectively means that from 
each magnetic map we are subtracting off any excess net flux, 
^ Br{9,<f>) sin 9 d9 d(j), a practice which leads to the introduc- 
tion of small errors in the resulting analysis. However, these 
errors are deemed to be less important than the inaccuracies 
resulting from the less-than-perfect knowledge of the radial 
magnetic flux on the Sun, including effects due to evolution 
and temporal sampling throughout each Carrington rotation 
and due to the lack of good radial field measurements of the 
flux in the polar regions of the Sun. 



^Tables of gy and fo™ are available from the WSO webpage at 
|http : //wso . Stanford ■ edu /Harmonic . rad/ ghli st . html 
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3. Multipolar Expansions and Their Evolution as a 
Function of Cycle 

3.1. Dipole Field (Modes with £ ^ 1) 

The solar dipolar magnetic field can be analyzed in terms 
of its axial an d equatorial harm onic components. As has long 
been known dHoeksemalllQSl . the axial dipole component, 
having a magnitude of |, is observed to be largest during 
solar minimum when there is a significant amount of magnetic 
flux located at high heliographic latitudes on the Sun. These 
two so-called polar caps possess opposite polarity, and match 
the polarity of the trailing flux within active regions located 
in the corresponding hemisphere that emerged during the pre- 
vious sunspot cycle. Long-term observations of surface-flux 
evolution indicate that a net residual amount of such trailing- 
polarity flux breaks off from decaying active regions and is 
released into the surrounding, mixed-polarity quiet-sun net- 
work. This flux is observed to continually evolve as flux ele- 
ments merge, fra gment, and move arou nd in response to con- 
vective motions (" Schrijver et al. 1997 1). but the long-term ef- 
fect is that the net residual flux is slowly advected poleward by 
surface meridional flows. Such poleward advection results in 
a net influx of trailing-polarity flux into the higher latitudes. 
At the same time, an equivalent amount of leading-polarity 
flux from each hemisphere cancels across the equator, as is 
necessary to balance the trailing-polarity flux advected pole- 
ward. Over the course of a sunspot cycle, this process is re- 
peated throughout subsequent sunspot cycles, during which 
flux from the trailing polarities of active regions eventually 
cancels out the polar-cap flux left over from previous cycles. 
Once the leftover flux has fully disappeared, the buildup of 
a new polar cap having the opposite polarity occurs by the 
subsequent activity minimum. 

In contrast to the axial dipole component, the equatorial 
dipole components, having magnitudes ^| and \Bl\, are 
largest during maximum activity intervals and weakest during 
activity minima. Individual active regions on the photosphere 
each contribute a small dipole moment that, aside from the 
small axial component arising from the Joy's Law tilt, is ori- 
ented in the equatorial plane. Together the equatorial dipole 
moments from the collection of active regions add vectorially 
to form the overall dipole moment. When many active regions 
are on the disk, it thus foUows that the equatorial dipole mode 
is likely to have a higher amplitude. During periods of quiet 
activity with few active regions on disk, the equatorial dipole 
ampUtude is minimal. 

Because the WSO data span three sunspot activity cycles, 
a bit of historical perspective on the evolution of the dipole 
can be gained, as shown in Figure [T] Figure [Tta) shows the 
amplitude of the axial dipole moment since mid 1976 and its 
evolution as a function of the activity level, represented in 
the figure by the sunspot numbej^ (SSN). It is also evident 



that, during the most recent minimum prior to Cycle 24, the 
magnitude of the axial dipole component is much lower than 
during any of the three previous minima (i.e., those preceding 
Cycles 21-23). The connection between the axial dipole com- 
ponent and the flux in the northern hemisphere is illustrated in 
the time-history of the mean flux density as integrated over the 
northern hemisphere, shown in Figure |2] We observe a slight 
lag for the northern hemisphere averaged magnetic flux with 
respect to the axial dipole coefficient due to the contribution 
of other axisymmetric modes possessing a different phase. 

Figure \V[h) illustrates the magnitude of the equatorial 
dipole since mid 1976. In step with the relatively lower num- 
ber of active regions during Cycle 23 when compared with the 
maxima for Cycles 21 and 22, the equatorial dipole strength 
is found to be lower during the most recent maximum than 
during the maxima corresponding to Cycles 21 or 22. 

Given the variation in sunspot cycle strengths throughout 
the past few centuries, we suspect that cycle-to-cycle varia- 
tions in the magnitudes of the axial and equatorial modes are 
not unusual. Proxies of the historical large-scale magnetic 
field, such as cosmic -ray induced var iations of isotopic abun - 
dances measured from ice-core data ( Steinhilber et aljr2012 ). 
also show such longer-term variation and thus seem to be con- 
sistent with this view. Interestingly, the range over which the 
variation in the ratio of the energies possessed by the equa- 
torial versus the axial dipole components is about the same 
for the three sunspot cycles observed by WSO, as shown in 
Figure [ijc). Longer-term measurements of this ratio unfor- 
tunately are not available due to the lack of a sequence of 
long-term magnetogram maps to which the harmonic decom- 
position analysis outlined in Section|2]can be applied. 

3.2. Reversals of the Dipole 

The process by which old polar caps are canceled out and 
replaced with new, opposite-polarity polar caps, as described 
in the previous section, manifests itself as a change in sign of 
the axial dipole amplitude throughout the course of a sunspot 
cycle. Such dipole reversals for the past three sunspot cycles 
are shown in Figure[3] where the latitude and longitude of the 
dipole axis are plotted with time. It is found that the dipole 
axis spends much of its time in the polar regions, and for only 
about 12-18 months during these cycles is it located equator- 
ward of ±45°. 

During these reversals, which occur at or near maximal ac- 
tivity intervals, the energy in the dipole never completely dis- 
appears. We find that the reduction in the energy (i?")^ in the 
the axial dipole is partially offset by an increase in the energy 
{B^^)^ + (Bl)"^ in the equatorial dipole. This results in a 
reduction of the total energy X)m(^i )^ dipolar modes 
only by about an order of magnitude from its axial-dipole- 
dominated value at solar minimum, as shown in FigureUa). 



■^Sunspot numbers with slightly different calibrations are available from 
various sources worldwide. In this article, we use the indices pro- 



vided by the Solar Influences Data Center at the Royal Observa- 
tory of Belgium, whose sunspot index data are available online at 
|http : / /www ■ side ■ be/ sunspot- index-graphics /sidc_graphics ■ php] 
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Fig. 1. — Evolution of the dipole (£ — 1) modes, as char- 
acterized by the (a) axial dipole coefficient B^, (b) equato- 
rial dipole magnitude '^J (B^^)'^ + {Biy, and (c) the ratio of 

their energies [iB~^)^ + {Bl)'^]/{B^)^ for the WSO (black) 
and MDI polar-corrected (magenta) data sets. Panel (d) shows 
the monthly smoothed sunspot number (SSN) from Solar In- 
fluences Data Center at the Royal Observatory of Belgium. 
The WSO data have not been corrected for known satura- 
tion effects that reduc e the reported values by a factor of 1 .8 
dSvalgaard et al.lll978h . 




1980 1990 2000 2010 
year 



Fig. 2. — Northern hemispheric mean flux density (magenta) 
and axial dipole coefficient from WSO (black), illustrating the 
connection between the axial dipole and the flux in each hemi- 
sphere. The downward spike in the mean hemispheric flux 
density occurring in 2001 is likely related to WSO sensitivity 
issues occurring during that time period, and may not be real. 
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Figure [3] indicates that, during a reversal when the axial 
dipolar component is weak, the axis of the equatorial dipo- 
lar component wanders in longitude. This seemingly aimless 
wandering occurs because the longitude of the dipole axis is 
primarily determined by an interplay amongst the strongest 
active regions on the photosphere at the time of observa- 
tion. As older active regions decay and newer active regions 
emerge onto the photosphere, the equatorial-dipolar axis re- 
sponds in kind. 

3.3. Quadrupole Field (Modes with £ = 2) 

The evolution of the energy contained in the quadrupolar 
{£ — 2) modes exhibit much more variation than the dipole, as 
shown in Figure |4] As with the equatorial dipole components, 
all of the quadrupolar modes have more power during greater 
activity intervals than during quieter periods, as illustrated in 
the evolution of the various quadrupolar modes plotted in Fig- 
ure |5] Furthermore, when large amounts of activity occur, it 
is possible for the total energy J2mi^2^)'^ quadrupo- 
lar modes to be greater than the energy J2mi-^T)'^ the 
dipolar modes at the photosphere. The ratio between these 
two groups of modes is shown in Figure EJc), from which 
it is evident that during each of the past three sunspot cy- 
cles there have been periods of time when the quadrupolar 
energy exceeded the dipolar energy by as much as a factor 
of 10. The corona, in turn, reflects the relative strength of 
a strong quadrupolar configuration of photospheric magnetic 
fields by creating complex sectors and possibly multiple cur- 
rent sheets that extend into the heliosphere. One example of 
such complex field geometry is suggested by the potential- 
field source-surface model of Figure |6] where a quadrupolar 
configuration having an axis of symmetry lying almost in the 
equatorial plane is seen to predominate. 

3.4. Octupole Field (Modes with £ = 3) 

As with the quadrupole, the octupolar modes contain more 
power during periods of high activity and less power during 
minimum conditions, as illustrated in Figure|2l The exception 
is the axial octupolar coefficient B^, plotted in panel (a) of 
Figure |7] which is nonzero during solar minima and exhibits 
sign reversals during sunspot maxima in a manner similar to 
the axial dipole coefficient B^. 

The behavior of the various m = modes can be under- 
stood by considering their functional symmetry: the func- 
tions are antisymmetric in 9 (i.e., antisymmetric across the 
equator) when the degree £ is odd, whereas for even £ the 
functions are symmetric in 9. The presence of polar caps 
during solar minimum, a highly antisymmetric configuration, 
is reflected in the similar evolution of the i?° and B^ coeffi- 
cients, which correspond to the axial dipole {£ = 1, m — 0) 
and octupole (£ = 3, m = 0) modes. The axial quadrupole 
{£ = 2, m = 0) mode does not share this behavior because, 
as a symmetric mode, it is not sensitive to the presence of the 
polar caps during solar minima. 




1980 1990 2000 2010 
year 



Fig. 3. — (a)-(c) MoUweide projections of the location of the 
dipole axis for the past three sunspot cycles (Cycles 21-23), 
as determined from WSO synoptic charts. The solid circles 
indicate the longitude and latitude of the dipole axis for each 
Carrington rotation, with every sixth Carrington rotation also 
indicated by an open circle. Grid lines (dashed) are placed 
every 45° in latitude and longitude for reference. The Car- 
rington longitudes of the central meridians of each projection 
are chosen to best illustrate the reversals, and differ in each of 
the panels. Panel (d) illustrates the latitude of the dipole axis 
as a function of time. The open circles in panel (d) correspond 
to same times as the open circles in panels (a)-(c). 
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Fig. 4. — Total energy in (a) dipolar modes X]m(^i")^' 
(b) quadrupolar modes J^mi^T)'^^ ^iid (c) their ratio 
T,miBT) / Em(^r)^ for the WSO (black) and MDI polar- 
corrected (magenta) data sets. Panel (d) shows the monthly 
smoothed SSN. The WSO data have not been corrected for 
known satura tion effects that reduc e the reported values by a 
factor of 1.8 dSvalgaard et al.lll978b . 




1980 1990 2000 2010 
year 

Fig. 5. — Evolution of the quadrupolar (£ = 2) modes, 
as characterized by the (a) axial quadrupole coefficient 
along with higher-order magnitudes of the (b) m — 



±1 modes 



and (c) TO — ±2 modes 



(5^2)2 + {Bj)^, for the WSO (black) and MDI polar- 
corrected (magenta) data sets. Panel (d) shows the monthly 
smoothed SSN. The WSO data have not been corrected for 
known satur ation effects that reduc e the reported values by a 
factor of 1.8 dSvalgaard et al.lll978h . 
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Fig. 6. — Representation of the coronal magnetic field in 
October 2000 for which the large-scale field is predomi- 
nantly quadrupolar This field is the result of a potential- 
field source-surface extrapolation (ISchatten et al.ll 1969) with 
an upper boundary of 2.5 Rq at which the coronal field is as- 
sumed purely radial. Both closed (black) and open (magenta 
and green, depending on polarity) field lines are shown in the 
model. Also shown is the contour of Br = at i? = 2.5Rq 
(thicker black line). 



The dependence of the coefficients on the degree £ is 
illustrated in Figure[8] where the time-averaged energies from 
the MDI data (spanning Solar Cycle 23) as a function of de- 
gree £ are plotted. Prior to averaging, the spectra were placed 
in two classes: Carrington rotations for which the SSN is rel- 
atively large (defined as when SSN>100) and rotations for 
which the SSN is relatively small (defined as when SSN<50), 
thus capturing the state of the Sun when it is either overtly 
active or overtly quiet. The figure indicates that the even- 
odd behavior is more pronounced during quiet periods, and 
these occur near and during solar minimum when the polar- 
cap field is significant. During active periods the even-odd 
trend is still recognizable, but because the polar caps are weak 
and the active regions are primarily oriented east-west (i.e., in 
the equatorial plane and thus contributing little power to the 
axisymmetric modes) the even-odd trend is less pronounced. 
We will further discuss the behavior of axisymmetric modes 
in the context of Babcock-Leighton dynamo models in Sec- 
tion 1431 



3.5. Full Spectra and Most Energetic Modes 

One property of the spherical harmonic functions IJ" (9, (p) 
is that the degree £ is equal to the number of node lines (i.e., 
contours in 9 and where ¥["■ = 0). In other words, the 
spatial scale represented by any harmonic mode (i.e., the dis- 
tance between neighboring node lines) is determined by its 
spherical harmonic degree £. As a result, the range of £ val- 
ues containing the greatest amount of energy indicates the 
dominant spatial scales of the magnetic field. 

To this end, we have averaged the non-axisymmetric power 
spectra from each of the datasets both over time and over m, 
and have displayed the result in Figure |9] As with Figure |8] 
we have divided the spectra into active and quiet classes de- 
pending on SSN. In the figure, it can be seen that the magnetic 
power spectra form a broad peak with a maximum degree oc- 
curring at £p^,^^ «25, corresponding to a size scale of about 
360°/^p,^,^ «15° in heliographic coordinates. Stated another 
way, this indicates that much of the magnetic energy can be 
found (not surprisingly) on the spatial scales of solar active 
regions or their decay products. 

Energy spectra determined from WSO charts (not shown) 
do not show the same broad peak at £p^,^^ w25 as found in the 
curves from the MDI-derived data shown in Figure |9] This 
is an effect of the lower spatial resolution of the WSO mag- 
netograph (which has 180" pixels, and is stepped by 90" in 
the east- west direction and 180" in the north-south direction 
when constructing a magnetogram) versus MDI (which has 
a plate scale of 2" per pixel in full-disk mode). The WSO 
magnetograph, as a result, does not adequately resolve modes 
higher than about £ = 15, creating aliasing effects even at 
moderate values of £ in the energy spectra. Accordingly, as 
longer time series of data from newer, higher-resolution mag- 
netograph instrumentation are assembled, the high-£ behavior 
of the energy spectra (such as those shown in Figs. |8] and |9|l 
may change due to better observations of finer scales of mag- 
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Fig. 7. — Evolution of the octupolar (£ = 3) modes, as 
characterized by the (a) axial octupole coefficient B^, along 
with higher-order magnitudes of the (b) to — ±1 modes 

{B^)"^, (c) TO = ±2 modes 



and (d) to = ±3 modes W [B^ ) 
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{Bl)\ 



-2^2 



{BlY, 



for the WSO 



(black) and MDI polar-corrected (magenta) data sets. Panel 
(e) shows the monthly smoothed SSN. The WSO data have 
not been corrected for known satura tion effects that reduc e 
the reported values by a factor of 1 .8 (ISvalgaard et al ] |l978h . 
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Fig. 8. — Time-averaged energies in the axisymmetric modes 
)^ as a function of £ for MDI original (black) and polar- 
corrected (magenta) data sets, for more active conditions 
(solid lines; defined as when SSN>100) and for quieter pe- 
riods (dashed lines; defined as when SSN<50). The interpo- 
la tion scheme used to correct the MDI polar flux is described 
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modes X)m>o(^™)^ ^ function of £ for MDI original 
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for quieter periods (dashed lines; defined as when SSN<50). 
The interpolat ion scheme used to correct the MDI polar flux 
is described in I Sun etaL (2011). 
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netic field. 

3.6. Primary and Secondary Families 

The projection of the solar surface magnetic fields onto 
spherical harmonic degrees allows us to delineate the main 
symmetries of the magnetic field. As noted in ^ the har- 
monic modes can be classified as either axisymmetric (to ~ 
0) or non-axisymmetric (to 7^ 0). Separately, the harmonic 
modes can be either antisymmetric (odd i + to) or symmetric 
(even i + m) with respect to the equator dOubbins & ZhangI 
1993h . Some authors refer to antisymmetric modes as dipo- 
lar and symmetric modes as quadrupolar (presumably be- 
cause the axial dipole and quadrupole modes usually possess 
the most power in their respective categories), while others 
synonymously assign th ese modes to either the primary and 



secondary family (e.g., iMcFadden et alj|1991l when charac- 
terizing the Earth's magnetic field geometry), respectively. 
In this article, we adopt the primary- and secondary-family 
nomenclature when describing the equatorial symmetry be- 
cause this avoids the confusion that may otherwise occur 
when, for example, it is realized that the equatorial dipole 
mode (i — l,m — 1) formally belongs to the "quadrupolar" 
family of modes (since ^ + to is even for this mode). 

One important result put forward by the geomagnetic com- 
munity is that the relative strengths of the primary and sec- 
ondary families are different during geomagnetic field re- 
versals and excursions. During reversals, the modes associ- 
ated with the secondary family predominate over primary- 
family modes, and during excursions this is not the case 
( Leonhardt & FabianI 12007 ). We now investigate whether 
analogous behavior is occurring in the solar setting, by deter- 
mining which harmonic modes are most correlated with the 
axial dipole and axial quadrupole. 

When applied to two variables, the Spearman rank correla- 
tion index p G [—1,1] indicates the degree to which two vari- 
ables are monotonically related. The index p is positive when 
both variables tend to increase and decrease at the same points 
in time. A rank correlation analysis is more general than 
a more common Pearson correlation analysis, which specif- 
ically measures how well two variables are linearly related, 
whereas the rank correlation analysis enables a determination 
of whether the time evolution of two mode amplitudes follow 
a similar pattern in time without regard to their (unknown) 
functional dependency. 

In Table [T] we list the degrees £ and orders to correspond- 
ing to the harmonic coefficients {t) that have the highest 
p (positive correlation) when compared with the axial dipole 
and axial quadrupole coefficients B'^{t) and ^"(i) (which 
peak at different phases of the sunspot cycle). The corre- 
sponding harmonic modes comprise the strongest modes in 
the primary and secondary families, respectively. We find 
that, among the mode amplitudes that are positively correlated 
with Bi, two out of three belong to the primary family. Simi- 
larly, for B2, 7 of the 10 most-correlated modes are members 



of the secondary family. 

These correlations indicate a preference in the solar dy- 
namo, at least as inferred from its surface characteristics, for 
modes belonging to the same family and thus having the same 
north-south symmetry characteristics to be excited nearly in 
phase. This preference is demonstrated further in Figures [TOl 
and[TTl in which the long-term trends of the time evolution 
of the first several axisymmetric mode coefficients are shown, 
after smoothing with a boxcar filter having a width of 1 yr. 
(We focus here on the axisymmetric mode properties because 
these modes are the only ones considered in most mean-field 
dynamo models, as discussed further in ^4.21 ) In Figure [TOl 
there is a clear correlation amongst the first few odd-£ and 
amongst the first couple of even-£ mode coefficients, a trend 
which is emphasized in Figure [TT] in which these same mode 
coefficients are overplotted. We note that the mode group- 
ings are not precisely in phase, as evidenced for example by 
the lag in ^ = 3 and especially the £ — 5 modes reversing 
signs with respect to the £ = 1 mode. When £ > 6 or so, 
these trends become much weaker amongst the axisymmetric 
modes (although Table [T] indicates that this is not necessarily 
true for the non-axisymmetric modes), presumably because as 
smaller and smallelr scales are considered the effects of the 
global organization associated with the 1 1-year sunspot cycle 
are less important in structuring the surface magnetic field. 

Figure [TT| additionally illustrates that the modes of the 
secondary family attain amplitudes of about 25% of the pri- 
mary mode amplitudes. Furthermore, the primary and sec- 
ondary mode families are out of phase: during reversals the 
primary modes become weak at the same time as the ampli- 
tudes of the modes associated with the secondary family be- 
come maximal, which was shown previously in Figure Uc). 
This same pattern is observed to occur during reversals of the 
axial dipole field of the geodynamo. As in the geodynamo 
case, we ascribe the relative amplitudes and phase relation 
between the primary and secondary families observed during 
solar dipole reversals as a strong indication that the interplay 
of the mode families play a key role in the process by which 
the axial dipole reverses. Hence, any realistic model of the 
solar dynamo must excite both families of modes to similar 
amplitude levels, and must exhibit similar coupling between 
modes belonging to the primary and secondary families. 

4. Theoretical Implications for Solar Dynamo 

As demonstrated in previous sections, the amplitudes of 
the various harmonic modes of the solar magnetic field are 
continually changing. During reversals, as the axial dipole 
necessarily undergoes a change in sign, other modes predom- 
inate such that the amplitude of the solar magnetic field never 
vanishes during a reversal. As a result of such reversals oc- 
curring in the middle of each 11 -yr sunspot cycle, during 
the rising phase of each cycle the polar field s and the emer 



gent poloidal fields h ave opposite polarity dBabcocklll961 



Benevolenskaval2004 ), while in the declining phase the polar- 
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Fig. 11. — Overplotted coefficients from Figure[TO]of the 
first 3 odd {£ = {1, 3, 5}; dark red, red, light red lines, respec- 
tively) and even (I — {2,4,6}; dark blue, blue, light blue, 
respectively) axisymmetric modes, as calculated from WSO 
synoptic maps. 
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Fig. 10. — (a)-(h) Coefficients for the axisymmetric 
modes of the first eight degrees ^ as a function of time, as cal- 
culated from the WSO (black) and MDl polar-corrected (ma- 
genta) synoptic maps and after boxcar smoothing with a width 
of 1 yr. Panel (i) shows the monthly smoothed SSN. The WSO 
data have not been corrected for known saturation effects that 



reduc e the reported values by a factor of 1.8 (ISvalgaard et al 
1978h . 
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Table 1: Ranking of the most positively correlated modes 
within the primary and secondary families for which £ < 9. 
The basis for comparison in each family is the lowest-degree 
axisymmetric mode belonging to each family, specifically the 
axial dipole (£ ^ 1, m — 0) and axial quadrupole (£ = 2, 
TO = 0) modes for the primary and secondary families, re- 
spectively. The most correlated mode is the next axisymmet- 
ric mode in each family. The equatorial dipole mode (£ = 1, 
TO = 1) is more correlated with the axisymmetric quadrupole, 
as expected from its symmetry properties. We note the pres- 
ence of 4 sectoral modes (for which ^ = m) in the list of the 
secondary family. 

ity of sunspots and active regions are aligned with the newly 
formed polar-cap field. 

We have already noted how the temporal modulation of 
the large-scale harmonic modes comprising the primary and 
secondary families during polarity reversals appears similar 



to that of the magnetic field of the Earth jMcFadden et al 
19911 iLeonhardt & Fabianll2007l) . We have illustrated that as 
the magnitude of the primary-family mode amplitudes (pri- 
marily those of the axisymmetric odd-£ modes _B^) lessen, 
the secondary-family mode amplitudes (particularly those of 
the equatorial dipole i?J and axial quadrupole B2) simulta- 
neously increase. Once the secondary-family modes have 
peaked, the primary-family modes grow as a result of the 
growing polar-caps. Such interplay between primary and sec- 
ondary families provides insight toward an understanding of 
the processes at play in the solar dynamo that are assumed to 
be r esponsible fo r the occurrence of the observed cyclic activ- 
ity (lTobiasll2002b . 

The presence of power in members of both the primary 
and secondary families indicates that the solar dynamo excites 
modes that are both symmetric and antisym metric with re- 
spect t o the equator. As was demonstrated by Robe rts & Stix 
(Il972h . this cannot occur unless nonlinearities exist or un- 
less basic ingredients of the solar dynamo (such as, for ex- 
ample the a and/or w effects, or the meridional flow) pos- 
sess some degree of north-south asymmetry. In light of the 
the parameter regime in which the solar dynamo is thought 



to operate, including large fluid and magnetic Reynolds num- 
bers Re and Rm believed to chara cterize the solar convec- 
tion z one (both of order lO^^-lO^^; IStixll2002l: lOssendrijver 



20031), one expects the Sun to possess a nonlinear dynamo. 
Detailed observations of the magnetic field in the solar inte- 
rior where dynamo activity is thought to occur are not avail- 
able, but the observed magnetic patterns and evolution pro- 
vide circumstantial evidence of turbulent, highly nonlinear 
processes that lead to complex l ocal and nonlocal cascades 



of energy and rnagnet i c helicity dA lexakis et al.ll2005L 12007 



Livermore et al. 2010t Strugareket al.. .2012 ). Yet, the pres 



ence of regular patterns formed by the emergent flux on the 
solar photosphere, as codified by Hale's Polarity Law, Joy's 
Law for active-region tilts, and the approximately regular cy- 
cle lengths, suggest that some ordering is indeed happening 
in the solar interior 

With the aim of distilling the necessary elements of the 
various nonlinear dynamos into a manageable framework, 
multiple authors have created idealize d models of the so - 

lar dynamo, including (for examplet Weiss et al ] (11984'); 

[pevnman & Gabriel (1990); Ruzmaikin et al.l(ll992l)~ [kiiobloc h& Lands 
( 1996|) ;lTobias ( 1997); Knoblochet al. ( 1998); M elbourne et all 
i200ll) ; 1Weiss & Tobias,~ 32000l) ; iSpiegeL (^Q^T Simi- 
larlV j for the g eodynamo t here a re many efforts , includ- 
ing 



Glatzmaier & R oberts (1995); Heimpel et al. (2005^1 



IChristensen & Aubert (2006); Busse & Simi tev.(,2008,);,Nishikawa & Ku 



([2008); Takahashi et al. (2008); Christense n^al.1 d2010l). A 



compl etely different approach h as been taken by Petrelis & Fauve 
( l2008h and lPetrelis et al.l(l2009b . who have developed simpli- 
fied models of the von Karman s odium (VKS) laboratory 
experiment (IMonchaux et al.ll2007h . In all of these idealized 
models, the modulations resulting from the coupling between 
magnetic modes from the differe nt families, or between the 
magnetic field and fluid motions (lTobiasll2002l) can be ana- 
lyzed in terms of the equations that describe the underlying 
systems. The variability of the most prominent cycle period 
develops as a result of the coupling of modes introducing a 
second time scale into the dynamo system, often leading to a 
quasi-periodic or chaotic behavior of the magnetic field, cycle 
length, and/or dominant parity. One can further understand 
via symm etry considerations how reversals and excursions 
arise (Gu bbins & Zhanglll993b . 



4.1. Reversals and Coupling between Modes 

To illustrate such a dyn amical system, following the work 
of IPetrelis & Fauvd (l2008h . we assume that the axial dipole 
and quadrupole modes are nonlinearly coupled. We can then 
define a variable A{t) — [t) + (t), where we have used 
the time-varying mode coefficients defined in Equation[T] and 
write an evolution equation that satisfies the symmetry invari- 
ance found in the induction equation, i.e., B — s> B. It then 
follows that the symmetry A — —A must also be satisfied, 
and that such a equation to leading order is 



dA 
dt 



/ii A + /i2 A + i/i + 1/2 A^A + 1/3 AA^ + , (5) 
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where fXi and i/i are complex coefficients and A = Bi — iB^ 
is the complex conjugate of A, and the quadratic terms 
have van ished due to symmetry con siderations. As dis - 
cussed in IPetrehs & Fauvei ('2008') and 'P etrehs et all (|2009|) . 
such dynamical systems are subject to bifurcations. In par- 
ticular, they demonstrate that this dynamical system can 
be characterized by a saddle-node bifurcation when com- 
paring its properties with so-ca lled normal form equations 
dGuckenheimer & Holmeslll982h . In such a bifurcated sys- 
tem, both stable and unstable equilibria (fixed points) exist, 
as illustrated in the left panel of Figure [12] For instance, 
if the solution lies at a stable point (for example, where the 
dipole axis is oriented northward), fluctuations in the system 
may disturb the equilibrium and push the magnetic axis away 
from its stable location. If such fluctuations are not strong 
enough, the evolution of the dynamical system resists the de- 
terministic evolution of the system and the system returns to 
its original configuration (in the example, resulting in an ex- 
cursion of the dipole), such as seen in the geomagnetic field 
dLeonhardt & FabianI 120071) . If instead the fluctuations are 
large enough to push the system past the unstable point, the 
magnetic field then evolves toward the opposite stable fixed 
point allowed by J5 — > —B (in the example, resulting in a 
reversal that changes the dipole axis to a southward orien- 
tation). Such behavior is also seen in the VKS experiment, 
from which is observed irregular magnetic activity combined 
with both excursions and reversals. Reversals result in an 
asymmetric temporal profile, with the dipole evolving slowly 
a way from its equilibr ium followed by a swift flip (cf.. Fig. 3 
of lPetrelis et aDl2009l) . 

In Figure[T3]we have overplotted the last three sunspot cy- 
cle reversals such that the zero-crossings of the axial dipole 
coefficients B^ for each cycle are aligned. It has been recently 
shown that the 10 major geomagnetic reversals for which 
detailed records exist occurring during the past 180 Myr 
possess a characteristic shape upon suitable normalization 
( Valet & FournieJ2012 ). This shape can be described as com- 
prising a precursory event lasting of order 2500 yr, a quick 
reversal not exceeding 1000 y r, and a rebound event of order 
2500 yr. IPetrelis et all (l2009l) show that the magnetic field in 
a simplified VKS laboratory experiment exhibits differing be- 
havior during reversals and excursions. During reversals, the 
magnetic field has an asymmetric profile that contains a slow 
decrease in the dipole, followed by a rapid change of polarity 
and buildup of the opposite polarity, whereas excursions are 
more symmetric. Additionally, after reversals the magnetic 
dipole overshoots its eventual value before settling down, 
whereas durin g excursions no such overshooting is measured 



(see Fig. 3 of IPetrelis et al. I I2OO9I) . For the solar cases dis- 



played in Figure [T3] we find that only the (green) curve of 
the reversal of cycle 22 exhibits an overshoot, whereas the 
other two cycles do not. Further, the rates at which the solar 
dynamo approaches and recovers from the reversal appear to 
be equal, leading to a symmetric profile, in contrast with the 
VKS results. Therefore, the Sun seems to reverse its mag- 



netic field in a less systematic way than other systems that 
have shown such behavior 

Analyzing such systems from a dynamical systems per- 
spective, when changing the control parameter (here, i?,„) 
past the bifurcation point, the stable and unstable points coa- 
lesce and merge and the saddle nodes disappear, as shown in 
the right panel of Figure [121 This act transforms the system 
from one containing fixed points to one containing limit cy- 
cles with no equilibria (e.g., Guckenheimer & Holmes 1982|), 
yielding an oscillatory solution that manifests itself as cyclic 
magnetic activity. Typically, large fluctuations are required 
in order to put the dynamical system above the saddle-node 
bifurcation threshold. 

In the case of the Sun, both the primary and secondary 
families are excited efficiently and a strong coupling betwee n 
them is exhibited. The model of IPetrelis & Fauve (I2OO8I) 
in spir it very close to the studies of Knob loch & Landsberj 
1996h or lMelbourne et al. I dlooi), may be used to guide our 



interpretation of the solar data. As illustrated in Figure[TT] the 
axisymmetric dipole and quadrupole are out of phase, such 
that their coupling may lead to global reversals of the solar 
poloidal field. To the best of our knowledge, however, the so- 
lar dynamo does not exhibit excursions of its magnetic field 
(unlike the geodynamo) but instead undergoes fairly regular 
reversals that take about one or two years to transpire (cf.. 
Fig. [3]). The solar dynamo is thus better approximated by a 
model in which a limit cycle is present. One may presume 
that the difference between the geodynamo and the solar dy- 
namo may be a result of the large degree of turbulence present 
in the solar convection zone, whereas the mantle of the Earth 
has a more laminar convective flow and thus is below the bi- 
furcation threshold where fixed points are still present. 

It may be the case that the solar dynamo is better described 
by a Hopf bifurcation, in which a limit cycle arises (branches 
from a fixed point) as the bifurcation parameter is changed. 
The dynamo instability that occurs as a result of the interac- 
tion of magnetic fields and fluid flows (such as auj d ynamos 



typically used to model the Sun, as summarized in iTobias 



1^ 

12002) often arises from a Hopf bifurcation. This allows the 
system to pass through domains having different properties, 
such as the aperiodic oscillations that characterize the grand 
minima and n onuniform sun spot cycle strengths of the solar 



dynamo (e.g., Spiege]||2009l and references therein). The data 



analysis shown here does not favor a particular type of bifur- 
cation, but does indicate efficient coupling between the pri- 
mary and secondary families. 

Yet another approach toward investigating magnetic re- 
versals is to develop detailed numerical simulations solv- 
ing the full set of MHD equations. Such three-dimensional 



numerical simulations in spherical ge 


iometry of the Earth's 


geodynamo dOlatzmaier & RobertsI 


19951 \ 


Lietal. 2002; 


Nishikawa & Kusanoll2008l: lOlson et < 


1. 12OII 


) or of the so- 


lar global dynamo ( Brun et alj 20041 


Browning et al. 20061 



Racine et al.ll201 ll) have looked at the behavior of the polar 



dipole vs multipolar modes. Even though such models have 
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large numerical resolution and thus possess a large number 
of modes, all have the property that the dominant polarity 
of the magnetic field follows the temporal evolution of a 
few low-order modes, even if in some cases the magnetic 
energy spectrum peaks at higher angular degree I. These 
findings suggest that the coupling between the primary and 
secondary family remains an important factor in character- 
izing polarity reversals for these simulations, and is thought 
to be Unked to a symmetry-breaking of the c onvective flow 
jNishikawa & Kusanol2008l: IOlson et al."201 1). Some studies 
of the geomagnetic field (e.g.. IClement.2004) even advocate 
for a coupling between two modes of the same primary fam- 
ily, such as the axial dipole and octupole. While in the solar 
data these modes are well correlated, the coupling between 
the primary and secondary families of modes seems more 
likely to be at the origin of the reversal, as demonstrated in 



Mean-Field Dynamo Models and the Axisymmetric 
Modes 



4.2. 



Mean-field dynamo models are found to capture the 



essence of the la r ge scale solar dynamo (iMoffat 



Ossendriiveii 120031: iDikpati et al.ll2004t ICharbonnea 



19781: 
20101) . 



At present, the most favored model is the mean-field Babco ck 



Leighton (BL) dynamo model (e.g.. IDikpati et al.l 120041) . in 
which the mean magnetic induction equation is solved us- 
ing empirical guidance for both the differential rotation and 
meridional circulation profiles, as well as for parameteriza- 
tions of the a-effect and poloidal-field source terms. In this 
section, we use the Stellar Elements (STELEM; T. Emonet & 
P. Charbonneau, 1998, priv. comm.) code (see Appendix A 



of iJouve & Brunll2007l for more details) to solve the axisym- 
metric BL dynamo equations, and investigate some of the 
consequences of the coupling between modes from the pri- 
mary and secondary families. In the interest of brevity, we 
refer interested readers to Appendix lA.ll for a listing of the 
governing equations associated with BL dynamos. 

In many BL solar dynamo models, the parameters govern- 
ing the imposed flows and the poloidal-field source terms are 
chosen based on their solar counterparts. When carefully cho- 
sen these terms favor modes from the primary family and thus 
antisymmetric field configurations, since this is what the Sun 
apparently favors much of the time. This is a result of the 
commonly used latitudinal profiles of the key dynamo ingre- 
dients (symmetric large-scale flows and antisymmetric alpha 
effect) combined with the parity in the BL mean-field dynamo 
equations, leading to a situation where modes of the primary 
family remain uncoupled to modes of the secondary family 
that aflows both dynamo families to coexist without much 
interaction. We consider the symm etry of the BL equati ons 
used here in Appendix | A .2l(see also iRoberts & Stixlll972 and 
Gubbins& Zhang! 1993 for broader discussions on this topic). 

To demonstrate these characteristics, we now consider a 
typical solution of the standard BL mean-field dynamo as cal- 
culated by STELEM. Figure[T4]presents the time evolution of 
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Fig. 12. — Schematic diagrams (adapted from ^ Petrelis et al 



E009) of a magnetic dynamo system on either side of a saddle- 
node bifurcation, with two distinct polarity configurations 
represented by B and —B. The coordinate axes represent 
states where the primary (as represented by the axial dipole 
U) or secondary (as represented by the axial quadrupole Q) 
families are dominant. In the left-hand diagram, stable (±i?s) 
and unstable (i-Bu) states present during the system's evo- 
lution are indicated by blue circles and red squares. Pertur- 
bations away from a stable point can either cause the system 
to evolve to the opposite stable configuration (if the perturba- 
tion is strong enough) or simply cause an excursion in which 
the system returns to the same stable state. In the right-hand 
diagram, corresponding to the same system at a higher mag- 
netic Reynolds number Rra the stable and unstable points 
have merged, and the system simply oscillates between the 
two configurations in a limit cycle. 
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Fig. 13. — The reversals, as defined by the magnitude of the 
axial dipole component for WSO, for the past sunspot Cy- 
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the resulting magnetic field patterns, and is thus analogous to 
the standard solar butterfly diagram. Performing a Legendre 
transform on the magnetic field reveals the degrees I of the 
dominant axisymmetric modes. Figure [15] illustrates that the 
odd i modes from the primary family dominate over the even 
ones by about five orders of magnitude in this model. This 
differs significantly from what is observed on the Sun, where 
the amplitude of the quadrupole is measured to be about 25% 
of the dipole amplitude for most of the time, becoming domi- 
nant only during reversals [cf., Fig.|4|c)]. The behavior of the 
standard BL model of Figure[T5]arises because of the symme- 
try characteristics of the BL dynamo equations. Because the 
model was initialized with a dipolar field, no modes from the 
secondary family are excited in the standard BL model shown 
in Figure [15] because no coupling exists between the primary 
and secondary families. 

If instead the calculation were initialized with a quadrupo- 
lar field (belonging to the secondary family), we find that the 
system eventually transitions to a state in which the primary- 
family modes predominate, as shown in Figure [16] The 
growth of the primary-family modes is due to the presence 
of a very weak dipole (likely of numerical origin) at the on- 
set of the simulations. In these models, the BL source term of 
Equation ( IA8I ) quenches the growth of the magnetic field once 
a certain threshold is passed, and as a result the maximum to- 
tal amplitude of the magnetic field is capped. The reason why 
the primary-family modes are preferred stems from the fact 
that the thresholds for dynamo action (based on the parame- 
ter Cs in Equation IIA4II ) are found empirically to be lower for 
the dipole than for the quadrupole {Cs ^ 6.12 vs. Cs ~ 6.25), 
meaning the dipole-like modes have a higher growth rate than 
the quadrupole-like modes. In this model, only briefly during 
the transition phase does the model possess a quadrupole of 
order 25% of the dipole, as in the Sun. 

Observations of solar photospheric fields, however, indi- 
cate that the Sun excites both families and does not strongly 
favor members of one family over the other, a situation that 
has apparently existed over many centuries. Even during the 
Maunder minimum, evidence suggests that this interval may 
have been dominated by a hemispherical dynamo with mag- 
netic activity located primarily in the southern hemisphere 
iRibes & Nesme-Ribes 1993), which can only be formed by 
a state in which primary- an d secondary-family modes pos- 
sess n early equal amplitudes ('Tobiad ll997l:lGallet & PetrelisI 
2009h . Consequently, the solutions presented in Figures [15] 



and [16] in which modes from only one family are preferred, 
are thus not a satisfactory model of the Sun. 

As advocated by Roberts & Stix ( 1972 ) and Gubbins & Zhan j 
d 1993b following their symmetry-based study of the solar 
dyn amo and the induction equation, and more recently by 
iNishikawa & Kus ano (2008) in their geodynamo simulations, 
a north-south asymmetry of the flow field specified in the BL 
dynamo, or alternatively an asymmetric poloidal-field source 
term, may allow the co-existence of both the primary and sec- 
ondary families. To investigate this effect we have performed 



two additional BL dynamo calculations, one with a BL source 
term and one with a meridional circulation that each generate 
both symmetric and antisymmetric fields (introduced via the 
parameter e of Equations ( |A9t and ( IA13I ) of Appendix I A. It , 
the combination of which yields asymmetric magnetic field 
patterns. 

We have run several dynamo cases, with the antisymme- 
try arising either in the BL source term or in the meridional 
flow profile, and with a range of amplitudes for the e param- 
eter from 10~^ to 10~^. All cases were initialized with a 
dipolar field. We find that when e is about 10^"^, the modes 
in secondary family grow until they reach about 35% of the 
dominant dipolar mode, as illustrated in Figures [TT] and [18] 
This result holds true regardless of whether the antisymmetry 
is introduced in the BL source term or in the meridional circu- 
lation profile, with very little difference in the resulting mode 
amplitudes. As expected, using a smaller e results in solu- 
tions where the primary-family modes dominate, while using 
a larger e yields a state where the secondary-family modes are 
comparable to the primary-family modes. Such results may 
indicate that Sun need only possess a weak degree of north- 
south asymmetry in order to behave as it does. 

5. Conclusions 

Cycles of magnetic activity in many astrophysical bodies, 
including the Sun, Earth, and other stars, are thought to be ex- 
cited by nonlinear interactions occurring in their interiors. Yet 
in some cases, such as the Sun, the cycles have approximately 
regular periods and in others, such as the Earth, there is no 
apparent periodicity. Dynamo theory indicates such a range 
of behaviors is expected, and whether the cycles are regular 
depends on magnetohydrodynamic parameters that character- 
ize the system, including fluid and magnetic Reynolds and 
Rayleigh numbers. As a consequence, the large-scale appear- 
ance of the magnetic field may provide clues toward the type 
of dynamo that may be operating. 

In this article, long-term measurements of the solar photo- 
spheric magnetic field are utilized to characterize the waves 
of dynamo activity that exist within the interior of the Sun. 
Synoptic maps from WSO (dating back to 1976) and MDI 
(spanning 1995-2010) are used to determine the spherical 
harmonic coefficients of the surface magnetic field for the past 
three sunspot cycles. We focus on the apparent interactions 
between various low-order modes throughout the past three 
sunspot cycles, and interpret these trends in the context of dy- 
namo theory. 

The multipolar expansions of the solar field as deduced 
from WSO and MDI data indicate that the axial and equatorial 
dipole modes are out of phase. During activity minima, the 
dipole component of the solar field is generally aligned with 
the axis of solar rotation, while the quadrupole component is 
much weaker During activity maxima, the dipole reverses its 
polarity with respect to the rotation axis, and throughout the 
reversal process there is more energy in quadrupolar modes 
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Fig. 14. — Latitude-time plots of (at the tachocline) and 
Br (at the surface) produced by a mean-field BL dynamo 
model that uses empirical guidance for the solar differential 
rotation and meridional flow profiles, and that is initialized 
with a dipolar magnetic field. The lower panel is analogous 
to the standard solar butterfly diagram. 




than in dipole modes. During the past three cycles, these re- 
versals have taken place over a time interval of about 2 yr 
to 3 yr on average. More indirect measures of solar activity, 
such as the sunspot number and proxies of the heliospheric 
field, seem to indicate that such regular activity cycles have 
persisted for at least hundreds of years with a period of ap- 
proximately 1 1 yr. The most recently completed solar cycle 
(Cycle 23) lasted for about 13 yr and while unusual, is not un- 
precedented. We note in passing that such modulations of the 
solar dynamo may be interpreted as a type of nonlinear inter- 
action betweenjhe turbulent alpha effect and the field and/or 
flows (lTobiasll2002h . 

The harmonic modes can also be grouped into primary and 
secondary families, a distinction that depends on the north- 
south symmetry of the various modes. For example, the ax- 
ial dipole harmonic is antisymmetric and is a member of the 
primary family. Alternatively, the equatorial dipole and ax- 
ial quadrupole modes are both symmetric with respect to the 
equator and thus are grouped together in the secondary family. 
When the evolution of the mode coefficients are analyzed in 
this way, we find that there is a trend for members of the same 
family to possess the same phasing, suggesting that modes in 
the same family of modes are either excited together and/or 
are more coupled when compared with modes of different 
families. This coupling is noticeable during reversals of the 
solar dipole, as less energy is present in primary-family modes 
than in secondary-family modes during these intervals. 

The historical record indicates that the geodynamo also un- 
dergoes reversals of its dipole axis (with respect to the rotation 
axis), but these reversals occur much more irregularly than in 
the solar case. Additionally, the dipole axis of the terrestrial 
magnetic field occasionally makes excursions away from the 
axis of rotation of the Earth, only to later return without ac- 
tually reversing. An examination of the large-scale harmonic 
modes of the geomagnetic field during these intervals indi- 
cates that the energy contained in secondary-family modes 
was significantly smaller during excursions than during re- 
versals. A strong quadrupole during geodynamo reversals is 
in line with the solar behavior; there is no parallel with excur- 
sions as excursions in the solar case have not been observed. 
Analogous behavior is observed to occur in the VKS labora- 
tory dynamo with respect to the relative strengths of the pri- 
mary and secondary families. 

We also examined the coupling of the mode families using 
a BL mean-field dynamo model computed using the STELEM 
code. Because of the symmetries in the magnetic induction 
and in the assumed profiles of the large-scale flow fields and 
BL source term, we find that the standard mean-field solar dy- 
namo model results in a state containing largely members of 
the primary family. This is a result of the dipole (a primary- 
family mode) being more unstable to dynamo action than the 
quadrupole. With a modest amount of asymmetry, imple- 
mented here either in the meridional flow profile or in the BL 
source term, we find from the models that both the primary 
and secondary families can coexist in the same model and in 
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the same proportions as in the solar dynamo. This can lead to 
a small lag between the northern and southern hemis pheres as 
is actually observed on the Sun dOikpati et al.ll2007 ). 
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Fig. 15. — Time evolution of the coefficients of the lowest- 
order harmonic functions of the surface magnetic field (as 
grouped by primary and secondary families) from the same 
BL dynamo model as shown in Fig. [14] In panel (a) are 
shown the evolution of the first several primary-family coef- 
ficients with ^ =1, 3, 5, and 7 in black, blue, green, and 
red, respectively. In panel (b) are shown the evolution of the 
secondary-family coefficients with £ =2, 4, 6, and 8, re- 
spectively in black, blue, green, and red. Panels (c) and (d) 
show zoomed-in sections of panels (a) and (b). 
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Fig. 16. — Time evolution of the same BL model as shown 
Fig.[T5](and using the same color scheme), but initialized with 
a quadrupolar magnetic field. 
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Fig. 17. — Time evolution of the same BL model as shown 
in Fig. [15] (and using the same color scheme), but with an 
asymmetric BL source term as implemented in Eq. iA9\ by 
setting e 10^^. 
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Fig. 18. — Time evolution of the same BL model as shown in 
Fig. [15] (and using the same color scheme), but with an asym- 
metric meridional flow profile as implemented in Eq. ( IA13I I 
by setting e = 10^'^. 
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A. Mean-Field Dynamo Formalism 



In this Appendix, we provide details regarding the Babcock-Leighton mean-field dynamo models discussed in 94.21 The pertinent 
equations are listed in ^A.ll and their symmetry properties are discussed in 9A.2I 

A.l. Mean-Field Equations 

Here, we briefly list the equat ions governi ng the axisymmetric mean-field dynamo models calculated in 94.21 A more detailed 
explanation can be found in, e.g.. |jouve et al.l (12008) . Following (.Moffatt.1978.) . the mean-field induction equation is 



dim 

dt 



V X {{V) X (B)) + V X {v' X b') -V X (ryV X (B)) , 



(Al) 



where the variables (B) and {V) refer to the mean parts of the magnetic and velocity fields, and v' and b' to their respective fluctu- 
ating components. The function rj is the magnetic diffusivity and is not necessarily a constant. The terms "mean" and "fluctuating" 
refer to the fact that a separation of scales has been performed, such that the mean quantities are computed by averaging over some 
appropriate intermediate size scale and the fluctuating quantities are the residuals. 

Working in spherical coordinates (r, 6, </>) and under the assumption of axisymmetry, we perform a poloidal-toroidal decomposi- 
tion and write the mean magnetic field B and mean velocity field V (for clarity the angle brackets ( and ) are omitted going forward) 
as 



B{r, e,t) = V X (A^e^) + 5^60 
V{r,9) = Vp + fir sin 60, 



(A2) 
(A3) 



where the poloidal streamfunction A^{r,9,t) and toroidal field Bci,{r,d,t) are used to generate B. The velocity field is time- 
independent, and is prescribed by profiles for the meridional circulation Vp{r, 9) and differential rotation f2(r, 9). 

Rewriting the mean induction equation (lAll i in terms of and B^, we arrive at two coupled partial differential equations for 

As and B^, 



dA^ 
dt 

dB^ 
dt 



v_ 




1 ^ 




Vt 








v_ 




1 ^ 




Vt 









R,^-V{wA^)^G,S 



1 d(wB^) d{j]lT]t) 



dr 



dr 



B,f, 
w 



(A4) 



ReB^V ■ Vp + Con7 [V X {A^e^)] ■ VVl, (A5) 



where w = r sin 9. The contribution to the transport term in the mean induction equation ( lAlb that arises from the fluctuating fields, 
namely the V X (v' x b') te rm, is p resent in the A^ equation above and in general is assumed to take a specific form in terms of 
the m ean magnetic field (cf., Babco ck 1961; Leighton 1969; Wang & Sheeley 1991; Dikpati & Charbonneau 1999; Jouve & Bru^ 



2007h . Here, we use a surface Babcock-Leighton (BL) term S{r, 6, B^) for this purpose which serves to generate new poloidal field. 



Additionally, equations (IA4l l and (IA5b have been nondimensionalized by using Rq as the characteristic length scale and Rq /rjt 
as the characteristic time scale, where rjt — 10^^ cm'^ s^^ is representative of the turbulent magnetic diffusivity in the convective 
zone. This rescaling leads to the appearance of three dimensionless control parameters Cji — iloRo/Vt^ Cs — soRQ/vt, and the 
Reynolds number R^ — uoi?©/?]*, where fio, Sf), and vq are respectively the rotation rate and the typical amplitude of the surface 
source term and of the meridional flow. 

Equations ( IA4b and iA5\ are solved with t he Stellar Elements (STELEM; T. Emonet & P. Charbonneau, 1998, priv. comm.) 
code (see Appendix A of Jouve & BrunI 2007 for more details) in an annular meridional plane with the colatitude e [0,7r] and 
the dimensionless radius r G [0.6, 1], i.e., from slightly below the tachocline (r w 0.7) up to the solar surface Rq. The STELEM 
code h as been thoroughl y tested and validated via an international mean field dynamo benchmarking process involving 8 different 
codes dJouve et al.ll2008h . At the latitudinal boundaries at 6* = and 6 — n, and at the lower radial boundary at r = 0.6, both 
Acf, and B^ vanish. At the upper radial boundary at r = 1, the solution is matched to an external potential field. Usual initial 
conditions involve setting a confined dipolar field configuration, i.e. Ac/, is set to (sin0)/r^ in the convective zone and to below 
the tachocline. To create the simulation shown in Figure [16] the simulation was initialized using a quadrupolar configuration with 
an Acf, of (3 cos 9 sin 9) j (2r"^) in the convection zone. In both cases, the toroidal field is initialized to everywhere. 

The rotation profile used in the series of models di scussed in this work ca ptures many aspects of the true solar angular velocity 
profile, such as deduced from helioseismic inversions ( Thompson et alj|2003 ). We thus assume solid-body rotation below r = 0.66 
and a differential rotation above this tachocline interface as given by the following rotation profile. 



nir,9) = n, + - 



1 +erf 



2(r - r,) 
di 



(ileq + a2 COS^ 9 + a4 CDs'* 9 — flc). 



(A6) 
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-0.136076 and = -0.145713. With this profile for n, the 



The parameters fieq = 1, = 0.93944, = 0.7, di = 0.05, 02 
radial shear is maximal at the tachocline. 

We assume that the diffusivity in the envelope 77 is dominated by its turbulent contribution, whereas in the stable interior rjc 't^ 
We smoothly match the two different constant values with an error function which enables us to continuously transition from r/c to 
Vt, 

iVt - Vc) 



ri{r) ^ric + 



1 + erf 



(A7) 



with 77c = 10^ cm^s-i and d = 0.03. 

In BL flux-transport dynamo models, the poloidal field owes its origin to the tilt of magnetic loops emerging at the solar surface. 
Thus, the source has to be confined to a thin layer just below the surface and since the process is fundamentally non-local, the source 
term depends on the v ariation of at the b ase of the convection zone. We use the following expression (which is a slightly modified 
version of that used in Jouve & Brunll2007h in order to better confine the activity belt to low latitudes: 



1 + erf 



r - r2 



1 - erf 



r - 1 



Bo 



COS0 sm^eB^{rc,0,t), (A8) 



where r2 = 0.95, o?2 = 0.01, Bq — lO^G. In the particular case of an imposed asymmetry between the north and southern 
hemisphere (Fig.fTTI) we introduce a modified source term, modulated by the parameter e, as follows. 



<S'asym('', 0, B^ 



1 + erf 



r -r2 



1 - erf 



B^{r,,e,t) 

Bn 



(cos 6* + e sin 61) sin^ eB^{rc,e,t). (A9) 



In BL flux-transport dynamo models, meridional circulation is used to link the two sources of the magnetic field, namely the 
base of the convection zone (where toroidal field is created via the latitudinal shear) and the solar surface (where poloidal field 
is introduced via the BL source term). In the series of models discussed in this paper, the meridional circulation is equatorially 
symmetric, having one large single cell per hemisphere. Flows are directed poleward at the surface and equatorward at depth (as in 
the Sun), vanishing at the bottom boundary at r = 0.6. The equatorward branch penetrates slightly beneath the tacho cline. To model 
the single cell meridional circulation we consider a stream function with the following expression ( Jouve et"aDl2008h . 



^(r,^) = -^^;-^sin 



7r(r - rb) 
1 - rb 



cos sin 6*, 



which gives, through the relation Vp = 'V X (V^e^), the following components of the meridional flow. 



(AlO) 



^(^-r»)ir-nr^.JAr-rby^ (3cos^^-l) 



nr (1 — VbY 

3r — Tfc . / 7r(r — rb) 
sm ' 



i-n 



i-n 

riT (r — Tb) 
" 1 - rfc (1 - n) 



7r(r - Tb) 

i-n 



2(1 - Tb) (r - n) 
nr (1 - rb) 



cos t) sine 



(All) 
(A12) 



with rb = 0.6. 

To introduce asymmetry into the model, an alternative to using the asymmetric source term of equation ( |A9t is to introduce an 
asymmetry into the meridional flow profile. Such an asymmetric meridional flow profile can be constructed using the following 
stream function. 



?/'asym(r, 9)^- ''^\ sin ( ] (cos + esm 0) sin 6, 



7r(l ~rb) \ 1 - rb 
which leads to the foHowing components of the meridional flow, 

2(1 - Tb) (r - rb)2 . fTi{r-rb) 



nr (1 — rb)' 



■ sm 



l-n 

again with = 0.6. 



3r — Tb . / 7r(r — rb) 
sm ' 



l-n 



rw {r — rb) 
1 - rfc (1 - rb) 



(3esin6lcos6l + 3cos^6l- 1) 



(A13) 



(A14) 



cos 



7r(r - rb) 

l-n 



2(1 - rb) (r - rb) 
Trr (1 - rb) 



(cos 6* + e sin 6*) sin 6*, (A15) 
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A.2. Symmetry Considerations 

Following iGubbins & Zhang ( 1993 ) it is straightforward to assess symmetry properties of various mathematical operators and 



equations. We adopt the superscripts ^ and ^ to indicate whether scalars or vectors are antisymmetric or symmetric across the 
equator, respectively. For example, products between a scalar and a vector of the form aF = G, where a and F are of like 
symmetry, yield a symmetric result (a^ F^' — > and a'^F"^ — !> G^'), whereas products between quantities of differing symmetries 
are antisymmetric {a'^F^ G^ and F^ G^). For the vector cross product F X G — H, when the two vectors F and 
G have the same symmetry properties the result will be antisymmetric (F^ X G^ — > and F^ X G^ — > H^), while the 
cross product between two vectors having opposing symmetries will yield a symmetric result (F^ X G^ H^). Additionally, 
the curl operator reverses symmetry (V X G^ and V X G^ — ?► H"^), while the Laplacian operator preserves symmetry 

With these properties established, the analysis of the symmetry properties of the magnetic induction equation, 

dB 

— X (V X B) + ijV^B, (A16) 

follows in a straightforward manner. For cases possessing a symmetric velocity field with respect to the equator, both terms on 
the right-hand side of equation ( IA16I ) preserve the symmetry of B. Thus, a dynamo having only a symmetric field B^ will remain 
symmetric over time, since both the transport term and the diffusion term of equation (IAI6I 1 generate symmetric field only. Likewise, 
a dynamo possessing an antisymmetric field B^ will preserve its antisymmetry over time. Because equation ( IA16b is linear in B, it 
follows that a magnetic field possessing mixed symmetry in the midst of a symmetric velocity field can be considered to be operating 
two independent, noninteracting dynamos: one that is symmetric and one that is antisymmetric. 

However, in cases with an antisymmetric velocity field V"^, the transport term on the right-hand side of equation ( IA16b provides a 
mechanism by which the symmetric and antisymmetric modes of B can couple. This coupling arises because an initially symmetric 
field B^ will first generate an antisymmetric field that in turn leads to the presence of both antisymmetric and symmetric fields, 
according to the right-hand side of equation ( IA16I ). Analogously, initializing with a purely antisymmetric field B^ wiU generate 
fields of mixed symmetry over time. 

This analysis procedure can further be applied to the mean-field induction equation. For example an analysis of equation (lAlb 
with an a-uj dynamo (i.e., where (^v' X 6') is set to a (B)), 

^ = V X {{V) x{B) + a {B)) r?V2 ((B)) (A17) 

S A 

indicates that, for an assumed symmetric mean velocity field {V) and an antisymmetric alpha effect a (which is the natural 
outcome of helical turbulence in a rotating fluid), such a rnean-fi eld dynamo will preserve th e symmetry (or antisymmetry) of 



the initial fields. Hence, as pointed out by Robe rts & Stix ( 1972 ) and McFadden et al. (Il99l l). a symmetric mean velocity field 



and an antisymmetric alpha effect do not couple magnetic field modes belonging to different families. Alternatively, if instead an 
antisymmetric mean flow (V)^ or a symmetric alpha effect are considered, this now enables a coupling between symmetric and 
antisymmetric mean fields. 

In a similar vein, the BL equations (IA4l l and iA5i . which are determined by performing the poloidal-toroidal decomposition on 
equation ( IA17b . can also be analyzed for symmetry. It is important to note that, by equation ( IA2I ). the poloidal streamf unction (^0) 
has a symmetry opposite to that of the mean magnetic field (B) it generates (and thus also to the corresponding toroidal field {B^)). 
We established above that the diffusion term in the mean-field equation ( IA17I ) preserves the symmetry of {B), and so it follows 
that the corresponding diffusion terms in the BL dynamo equations (IA4b and ( IA5b will serve to preserve the symmetries and 
B(f,. Likewise, because the large-scale transport term V X {{V) X {B)) term in equation ( IA17b preserves the symmetry of {B) 
whenever (V) is symmetric, it follows that the analogous terms in the equations (IA4b and (IA5b also preserve the symmetry of the 
system as long as the imposed velocity field is symmetric. For the BL dynamo considered here, an antisymmetric poloidal velocity 
streamfunction ^j^, as in equation ( lAlOb . yields a symmetric meridional flow profile Vp^ , since = V X (V'e,/,), which in turn 
gives a symmetric mean velocity {V)^ from equation (IA3b . Therefore, the imposed velocity field as defined by equations ( IA3b 
and (lAlOb will preserve the symmetries of and B^. Lastly, the source term S as defined by equation (IA8b also preserves the 
symmetry of A^, since it is comprised of a series of symmetric coefficients multiplied by cos9 B^. Therefore, an antisymmetric 
toroidal field implies a symmetric source term that in turn serves to preserve the symmetry of A^ (and thus (B)), and the same is 
true when the toroidal field is symmetric. 

For these reasons, the dynamo whose characteristics are illustrated in Figure[T5] which was initialized with a dipolar field (which is 
antisymmetric), preserves its antisymmetry with time since all of the terms in equations (IA4b and (IA5b preserve the initial symmetry. 
Indeed, the amplitude of the secondary-family modes (which are symmetric) remain low in this model, as shown in Figures \T5\ h) 
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and (d). In the dynamos whose characteristics are displayed in Fi gures [TTjorfTSl this effect is responsible for the growth of symmetric 
mean fields, even though both models were initialized with the same antisymmetric mean magnetic field. 

It is therefore a direct outcome of symmetry considerations that in standard mean-field dynamo models either one or the other 
families of magnetic fields is excited. In the experiments discussed earlier in 34.21 we controlled the degree to which the symmetries 
were mixed via the parameter e in equation iA9i and in equations (IA14l i and (|A15l l. which led to the dynamos illustrated in FiguresfTTl 
orfTSl respectively. In both cases, e was chosen to yield a dynamo where the end state contained secondary-family amplitudes of 
about 25%, as is observed on the Sun; other choices of e will lead to end states with different ratios. 
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